r 0 = 112 MeV. This is different than the value which 
would be found by taking the half-width from either of the 
plots in Figs. 1 or 2. Taking the half-width from the plot of 
Fig. 2 would be equivalent to finding the difference in cen¬ 
ter of mass energies at phase angles of 45° and 135°. With 
this data the result of doing this would be about 122 MeV 
(see Perkins 13 for an Argand plot which is used to find the 
half-width by this latter method). 

The value of T 0 could also be found from the slope of a 
plot of 8 (E ) versus center of mass energy in Fig. 3. Using 
L’Hopital’s rule on Eq. (11), one finds that r o is given by 
twice the inverse of the slope of Fig. 3 evaluated at E 0 and 
with 8 (E ) converted to radians. 

IV. DISCUSSION 

The delta 4 " + resonance in pion-proton scattering is a 
pedagogically good example: It is necessary to use relativis¬ 
tic kinematics in calculations, the approximate nature of 
the Lorentzian shape is illustrated, and it provides a clear 
example of how resonance parameters are determined by 
particle physicists. In addition, the important concept of 
isospin conservation can be illustrated by comparing the 
n + +p—* A + + cross sections to those for the 
tt~ -I- p—►A 0 . 14 If isospin conservation holds, V a for this 
latter reaction in Eq. (7) should be equal to T/3, so that the 
cross section as given by Eq. (10) would be a factor of three 
smaller for o\ir~ + p—* A 0 ). Carter et al. 15 give data for the 
(ir~ + p—* A 0 ) reaction as well as for the (n ' + p—+A + + ) 
reaction. Additional factors must be taken into account in 


this decay, however 15 : A small fraction of the decays take 
place by A 0 — >tz + y and the density of final states factor is 
slightly different for the decay A 0 —>7r _ + p than for the 
decay A 0 —>-77- + «. When these factors are included, a simi¬ 
lar analysis to find E 0 and r 0 for the delta 0 can be made. 
The importance of good data and its proper analysis would 
be seen in attempting to compare the rest mass energies of 
the delta + + and delta 0 . 
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Two types of damped harmonic motion are considered: motion that has damping terms 
proportional to the speed and the square of the speed and motion with damping that is constant 
and proportional to the square of the speed. These equations are solved analytically within the 
averaging approximation for small angle oscillations. Unknown coefficients in the equations of 
motion are then found by comparing their solutions with experimental data. Finally, the derived 
coefficients are used numerically to find values for 9 n0 , the angular speed of the pendulum at the 
lowest point of its (n + l)th swing, for a large release angle. The numerical solutions are compared 
with experimental results to validate the derived coefficients in the equation of motion. 


I. INTRODUCTION 

Damped harmonic oscillation is a source of fascination 
to many physics teachers. The equations of motion are 
slightly complex, and yet their solution is within the grasp 
of the average physics major. 

There are a number of experiments which mechanically 
illustrate damped harmonic motion. Most of these experi¬ 
ments involve the use of a pendulum; some measure the 
change of amplitude 1 and others the variation of the peri¬ 
od. 2 In the experiments described in this paper, we observe 
the decay of the speed of the pendulum at the lowest point 
of its swing. 3 This approach has only become practical with 


the advent of the microcomputer as a laboratory tool. Only 
a computer can react quickly enough to measure these 
successive speeds as the pendulum motion decays. In this 
paper we have considered pendula consisting of a rigid alu¬ 
minum bar supported on cone bearings, a golf ball, and a 
large polystyrene ball; both of the balls are supported by 
cotton thread. 

The computer used in this experiment is the Atari 800; 
however, the method described in Sec. II is quite general 
and can be easily implemented on many other microcom¬ 
puters. 4 ' 5 

In Sec. Ill we consider the equations of motion of a 
damped pendulum with small angle of swing and show 
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how the damping coefficients may be found from the exper¬ 
iment of Sec. II. Detailed solutions of the equations of mo¬ 
tion using the averaging approximation can be found in the 
Appendix. Sample experimental results are included and 
the associated damping coefficients found. These coeffi¬ 
cients are tested with data for large angles of the pendulum 
swing by comparing numerical solutions to the equations 
of motion that are predicted. The numerical method used is 
described in Sec. IV. 

II. EXPERIMENTAL METHOD 

A popular method of analyzing pendulum motion with 
the aid of a computer uses a potentiometer as the pendulum 
bearing. 6 As the pendulum swings, the varying resistance 
of the potentiometer is recorded by the computer. This 
method produces interesting results, but has the disadvan¬ 
tage of imposing a constant damping due to friction within 
the potentiometer. We found it possible to study a wider 
range of damping modes by measuring the speed of the 
pendulum at the bottom of its (n + l)th swing (<9„ 0 ) and 
observing how the ratio (0 nO /<9 (X) ) declines with the number 
of swings («). 

(9 n o) was determined from the time the pendulum took 
to pass through a photogate. This time may be measured in 
a variety of ways. We used a simple photogate interfaced 
through the game port to an Atari 800 microcomputer. 7 A 
short machine language routine 8 in the computer deter¬ 
mined the ratio \d nQ /0 m \ from the time that the photogate 
was blocked. The same routine was used in conjunction 
with the system clock to measure the small angle period of 
the pendulum. 

The photogate consists of an FPT 100 phototransistor 
and a 1-mW laser. The laser has two important advantages 
over conventional light sources. The laser may be set at 
some distance from the phototransistor, thus giving plenty 
of room for the pendulum and second, the tight beam of the 
laser obviates the need for collimation. 

Two of our pendula were a golf ball and a large polysty¬ 
rene ball, both of which were supported by cotton thread 
about 1 m in length. The thread was fastened to the golf ball 
with the aid of a small nail and to the polystyrene ball by a 
paper clip pushed into the polystyrene. The polystyrene 
ball had a diameter of 14.5 cm and was purchased from the 
craft supply section of a local hardware store. The third 
pendulum was a meter of 1/4-in-diam aluminum support¬ 
ed on cone bearings. 

The polystyrene ball was so large that \9 n0 /9 00 \ for the 
pendulum formed from it could not be found accurately by 
measuring the time the diameter of the ball took to pass 
through the photogate. This difficulty was resolved by 
passing the string supporting the polystyrene ball through 
a small metal cylinder of about 0.5-cm diam, which was 
allowed to rest on the ball. The metal cylinder was then 
used to cut the photogate rather than the diameter of the 
polystyrene ball. This solution was allowed because it was 
the ratio of angular velocities that was being measured. 

Once the apparatus was set up and the Atari suitably 
programmed, the experiment consisted of releasing the 
pendulum from a known angle and recording the various 
values of|0„ o /0J as the motion decayed. The period asso¬ 
ciated with each swing was also recorded. 

III. ANALYSIS AND RESULTS 


form for the equations of motion of various pendula having 
different damping terms. Initially, it was hoped that the 
large polystyrene ball would have damping proportional to 
the square of its velocity, 9 while the rigid pendulum would 
have a constant damping. The experiments bore out the 
first supposition, but found that the rigid pendulum also 
had a significant damping proportional to the square of the 
velocity. 

A general form for the equation of motion was assumed 
and solved using the approximation technique for small 
angles of swing. The solution was compared with the ex¬ 
perimental data gathered at small swing angles to see if it 
was of the predicted form and to determine numerical val¬ 
ues for the parameters of the equation of motion. Finally, 
the pendulum was released from a large angle and the ex¬ 
perimental data compared with data generated numerical¬ 
ly from the equation of motion. This last comparison was a 
sensitive test of the parameters found in the first part of the 
experiment. 


A. The ball pendula 


It is assumed that the equation of motion is of the form 

0 = -col sm9-K\9\9-J9 (3.1) 

which becomes, for small angle swings, 

e~-(ol 9-K\e\b-je. (3.2) 

In order to solve (3.1), we assume that the damping terms 
are small enough that the averaging approximation may be 
applied. The solution, derived in the Appendix [see (A13)] 

Hf 1 )) + (S 


X 



- 1 


(3.3) 


If the damping term J is small and we take a time t close 
enough to the beginning of the oscillation, we may ignore 
terms in Jt 2 and higher powers. Expanding the exponen¬ 
tials in (3.3) we find 

0~c o { 1 -(- [ J /I + 4A'w 0 a 0 /(37r)]t} -1 cos(a> 0 /). (3.4) 

This approximate solution should be valid for short 
times after the pendulum begins swinging and is justified by 
comparison with experiment. 

From (3.3) it is clear that when 9 — 0, t = (2n + l)7r/ 
(2<y 0 ), with n an integer. Thus the angular speed at the posi¬ 
tion 9 = 0 is given by 

9 n0 ~( - l)"KfloH 1 + [ J/2 + 4ATc[> 0 a 0 /(3jr)] 


X[(2n + \)ir/2co 0 )]} 1 . (3.5) 

Note that this speed is not necessarily the maximum speed 
of the pendulum and that the approximation will break 
down if n becomes too large. 

In the experiment we measure the ratio of the magnitude 
of the speed of the (n + 1 )th swing 1 9 n0 | to the magnitude of 
the speed of the first swing |£?oo|, with the speeds taken at 
the bottom of the swing, where 6 = 0. 

From (3.5) we can see that 

l^oo/^ol^l + kirn\co 0 [\ + (kir/.2a> 0 )]\~' , 
where 


The overall objective of the experiments was to find the 


k = J /2 + 4X<y 0 a 0 /(3flj , 


(3.6) 
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INITIAL ANGLE (DEG.) 



0 5 10 15 20 25 30 

SWING NUMBER (N) 


Fig. 3. Data for the large polystyrene ball released from an initial angle of 
70°. • are experimental points, while the + represent values derived 
numerically from (3.9). It is assumed in this case that the air resistance is 
only proportional to the square of the speed. 


Fig. 1. Data for the large (14.5-cm diam) polystyrene ball pendulum show¬ 
ing how.m, the slope of \6 M /8 m \ with swing number, depends on the 
angle of initial release. 

which is a straight line with a gradient m given by 
m = kv\co 0 [\ +kv/(2co 0 )] j -1 . (3.7) 

Expanding k and rearranging we find 
m /{2 - m) = (2K/3)a 0 + Jir/{Aa 0 ) . (3.8) 

t . . 

m, the gradient of the graph of 1 9 m /9 n0 1, and the number 
of swings are easily measured using the methods described 
in Sec. II. a 0 , the angle of release of the ball, is also mea¬ 
sured easily; thus Eq. (3 8) offers a practical method of de¬ 
termining the damping coefficients J and K. 

Plots of m/(2 — m) against a 0 for the polystyrene ball 
and the golf ball are shown in Figs. 1 and 2, respectively. 


From these graphs, values of K and J were determined for 
each pendulum by using Eq. (3.8). The following equations 
of motion for large angle swings were predicted: 

For the polystyrene ball, 

6= -og sin 6 -0.240 9\9\-, (3.9) 

for the golf ball, 

9 = — col sin 0 — 4.75 X 10~ 3 9 

- \.\9X\0~ 2 '9\9\ . (3.10) 

Numerical solutions generated from Eqs. (3.9) and (3.10) 
are compared with experiment in Figs. 3 and 4. The fit of 
the theory to experiment in both cases is quite good and 
justifies the approximations we made in setting up Eq. (3.8). 

It is interesting to note that the air resistance of the large 
polystyrene ball may be considered to be proportional only 
to the square of its speed. Derivation of this form of air 
resistance is a simple introductory mechanics exercise. 10 



O 10 20 

INITIAL ANGLE (DEG) 

Fig. 2. Data for a golf ball pendulum showing how m, the slope of 
I ^o/^ool with swing number, depends on the angle of initial release. 


B. The rigid pendulum 

The rigid pendulum is a meter length of aluminum rod 
pivoted on cone bearings. The cone bearings provide a con- 
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Fig. 4. Data for the golf ball released from an initial angle of 80°. • are 
experimental points, while the + represent values derived numerically 
from (3.10). In this case, the air resistance is assumed to be both directly 
proportional and proportional to the square of the speed. 


993 


Am. J. Phys., Vol. 53, No. 10, October 1985 


M. F. Mclnemey 


993 






stant damping term; however, this one damping term alone 
does not properly describe the motion. A damping term 
proportional to the square of the velocity must also be in¬ 
cluded. The equation of motion is thus of the form 

8 = -col sin e- k\ e\e- lb/\ e \. (3.11) 

If the damping constants are not too large and 0~sin 8, 
Eq. (3.11) may also be solved by the averaging method. The 
solution to Eq. (3.11), derived in the Appendix [see (A. 13)] 
is 

8 = cOq '(3 L /2 K) 112 cosK 1 ) tan 

X {tan~'[a 0 co 0 (2K /7>L ) 1/2 ] -2tn~\2LK/1) m \ . 

(3.12) 

From (3.12) we find 

\8 n o/8oo\ = tan | tan “ 1 [atfo 0 (2K /3Z. ) 1/2 ] 

- (2m + IK" \1LK/l) m ) 

X(tan(tan _1 [ao^o(2^/3Z, ) 1/2 ] 
-coo X (2LKn) U2 })~ 1 . (3.13) 

We now need to compare the expressions for j0 nO /0ool l n 
Eq. (3.13) with experiment to determine values of the 
damping constants J and K. The obvious method of finding 
/and K is to find the number of swings n 0 at which 8^/8^ 
= 0. In this case we see from (3.13) that 

tan - l [aco 0 (2K/lL ) 1/2 ] = (2n 0 + IK" '(2LK/l) in . 

(3.14) 

If (2n 0 + 1)K) _1 (2LA/3) l/2 is sufficiently small, we 
may expand the tangent to second order and find 

a 0 (2« 0 + \)~ 1 chLo)q 1 \ \ + (2 n 0 + l) 2 2LK/(9co 2 0 )] . 

(3.15) 

Unfortunately, we were not able to use this method in 
practice because the exact value of n 0 is difficult to mea¬ 
sure. Instead we note if <5 <jc, then 

tan(x + 35)/tan(x + 5)~ 1+25(1 + tan 2 *)/tan x . 

(3.16) 

This approximation may be applied to Eq. (3.13) taking 
n = 1: 


x^Um-'la^TK/lL) 112 ] 

and 

8={2LK/3) U2 /o ) 0 , 
to give 

8 io 2L _ 4Ka 0 

^oo (flo^o) 3 

The assumption underlying the approximation (3.16) 
may be expected to hold true in most cases of light damp¬ 
ing, where L and K are less than 1 and the constant friction 
term L is greater than the square term K. 

Experimental results are shown in Figs. 5 and 6. Values 
of \e i0 /e M \ were found for a number of initial release an¬ 
gles a Q and are plotted in Fig. 5. Also plotted in Fig. 5 is the 
least-square fit of Eq. (3.17) to the data, which gives the 
values K = 0.021 and L = 0.032. 

Using these derived values of K and/, the expected vari¬ 
ation of 1 8 n0 /8 a o | and the period with increasing number of 
swings from an initial release angle of 90° were found by the 
numerical solution 11 of Eq. (3.11). These results are com¬ 
pared with experimental data in Fig. 6. 

IV. NUMERICAL SOLUTION 12 

In the work described here, we used the fourth Euler 
approximation, 13 iterated as 

8 m = -co 2 0 sin 8 m —K \8 m \8 m -J8 m -L (8 m /\8 m \), 

8 m + i^8 m +8 m St, 8 m + , c^8 m + 8 m + x 8t, 

t m+l ~t m +8t, (4.1) 

with initial conditions t 0 = 0 ,8 0 — 0, and 8 0 = a 0 . We were 

interested in finding how \8„ 0 /8 m) \ and the period changed 
as the number of pendulum swings increased. To find both 
of these relations accurately using the Euler iteration 
would have required a very small time increment 8t, and 
consequently, several hours of microcomputer time. This 
inefficiency was avoided by taking advantage of the expect¬ 
ed nature of the solution. 

If the solution is a modified sine curve, which experi¬ 
ment suggests, then within a small enough time interval of 
8 = 0, the solution will be a straight line. In this case the 
time at which 8 = 0, t [8 = 0) may be derived by interpolat- 



Fig. 5. Data for the rigid pendulum from 
which the coefficients associated with 
the equation of motion (3.11) may be 
found. • are experimental points; the 
line is the least-square fit of the data to 
Eq. (3.17). 
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Fig. 6. Data from the rigid pendulum released from 90°. • are experimen¬ 
tal points, while + arederivedfromthenumericalsolutionof(3.17)using 
a time interval of 0.05 s. The period is the lower set of points. 


ing between the points 9 m and 9 m + l , which straddle 0 = 0. 
With this assumption 

t(0 = o)~t m +8t\e m /d m+ ,\{\ + \e m /9 m+i \)-' 

=t m + DA . (4.2) 

The computer found 0 m and 9 m + , by looking for 
9m*8 m + 1 <0- 

Interpolation was also used to determine the speed at 
6 = 0 and hence |0 reO /0oo|, although approximating the so¬ 
lution by a straight line is inappropriate in this case, and a 
fourth-order polynomial was used instead. 14 0„ o was found 
by using DA from Eq. (4.2) as the time argument of this 
polynomial. Finally, this careful calculation of t [6 = 0) 
and 0„ 0 allows us to apply a corrector to the Euler iteration 
at every half-swing by resetting the initial condition in Eq. 
(4.1) to t 0 = t (6 = 0), 0 O = 0„o, and 6 0 = 0. 

The method described here was tested on the undamped 
pendulum with a period of 1.4 s and St = 0.05 s, giving a 
numerical solution for 0„ o which was within 0.05% of the 
analytical value after 70 swings. 

V. CONCLUSION 

The experiments described in this paper offer students 
examples of real life damped harmonic oscillators that can 
be analyzed quantitatively. The experiments also introduce 
the student to simple numerical analysis and judicious use 
of approximation. The comparison of the observed decline 
of 9 n0 /9 00 with that predicted numerically is a satisfying 
way to see how “correct” the parameters are. 

We have considered motion with damping that is con¬ 
stant and proportional to the velocity and motion that is 
proportional to both the velocity and square of the velocity. 
A pendulum oscillating in a viscous liquid would have a 
damping which is the sum of constant terms and terms 
proportional to the velocity. We have not considered this 
case, but it is dealt with well by Ricchiuto and Tozzi. 1 

At Rose-Hulman, we are introducing the golf ball sec¬ 
tion of these experiments as part of our advanced laborato¬ 


ry. We hope that interested students will expand on the 
experiment, perhaps in the direction of a comparison of the 
air resistance at low speeds of different brands of golf ball. 

In general, the experiments can be extended to deter¬ 
mine the difference in air resistance of different shapes and 
the effect of air temperature and pressure on air resistance. 
There is scope for a number of student projects. 
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APPENDIX: THE AVERAGING METHOD 14 

The averaging method is a means of solving differential 
equations of the form 

6= -co\9 + /(0,0), (Al) 

where the damping terms /(0,0) are small, with little effect 
on the motion during any one period. 

If this assumption is true, then there is a solution of the 
form 

6 = a(t ) cos[u > 0 1 + ^(t)] , (A2) 

with 

J ~27r/co Q 

dt f(6,6)sm[o) 0 t+ <!>] (A3) 

o 

and 

if>(t )~ - (2tt-)- 1 " dt(f(9,9 ) cos . (A4) 

In working out these integrals, one makes the approxi¬ 
mation within the integral that a is independent of time. 
This approximation is simply a repetition of the assump¬ 
tion from which we began. 

The expressions (A3) and (A4) are derived in Baierlein. 


A. Ball pendula 

We will use the averaging method to solve the small an¬ 
gle equations of motion (3.2) and (3.11) beginning with (3.2), 
for which the equation for the ball pendula is 

9= -col 0-K\9\9- J9. (A5) 

We see that f[6,9 ) in this case is — (K\9 \9 + J9 ). From 
(A4) we find 

n2 v/(o 0 

— <b = (2ir)~ l dt col Ka\sm(a) 0 1 + <f>)\ 

Jo 


X sin(u > 0 1 + <b) cos(*u 0 1 + 

t'2ir/(o 0 


+ (277-) 


f 


dt Ja 0 


X sin(<u 0 1 + </>) cos (« 0 1 + <f>). (A6) 

Both terms in this equation are zero, as can be seen by 
substituting x = co 0 t + d>, whereupon (A6) becomes 

/•27T+ tf> 


- <t> = (277-) 


i 


dx u> 0 Ka|sinx|sinx cosx 


■c 


•2v + tf> 


+ (277-)“ 

= (2ir)-'Kaco () ^ 


dx J sin x cos x 
dx sin 2 x cos x 
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r*2ir /»2ir+ <j> \ / \ r2w/co 0 

j dx sin 2 x cos x + J dx sin 2 x cos xj 4 = — —j J dt \aa>l |sin(to 0 1 + 


+ (277-) 


"'(C 


dx sin x cos x 


I 


+ dx sin x cos x + 


I 


2v 4- <b 


dx sin x cos x 


(A7) 

If the integrals in (A7) between 2-ir and Itt + <f> are replaced 
by integrals between 0 and <ji it is easy to see that 4 = 0. For 
convenience we define the initial conditions such that 
4 = 0. 

From (A3) and by taking <f> = 0 we find that 

/*2tr/a> 0 


J ~zir/co 0 

dt Ka 2 col I sin co 0 t 
o 


sin 2 co 0 1 


+ 


J ~2ir/eo 0 

< 

0 


dt Jco<fl sin 2 co 0 t. 


(A8) 


In working out these integrals we note that 

J ~ 2 tt/co 0 rn/(o 0 

dt |sin(w 0 t)|sin 2 (<a 0 t) = 2 dt sin 3 (ty 0 f) (A9) 
o do 


and remember that in the averaging approximation, a{t ) 
may be considered a constant within the range of the inte¬ 
gration to obtain 

a(t) = - (4tf/3fl-Ktf 2 (t) - (J/2)a(t ). (A10) 


a(t) is found from (A10) by rearranging (A10) into the inte¬ 
gral form 


i>[& 2+ (f)r=-j> ,Ani 


The initial conditions are that at t = 0, a(t ) = a 0 . Inte¬ 
grating (A11) we find 

t = 2J~ x lnjaj J/2 + 4Kco<fl[t )/(3 v)] 


Xa^-'lJn + lKco^M]- 1 ] , (A12) 

which may be rearranged to give a(t ) and hence, using (A2), 
a solution for 6 : 


8 = (f) c “ K,, (( 7 )“ p (f) 

+(f)-Mf)-r' ,ai31 


B. Rigid pendulum 

The small angle equation of motion for the rigid pendu¬ 
lum is assumed to be (3.11) 

0= - col 6 ~ K\d\6 - Ld/\6 \ , (A 14) 

for which case/(0,0) is (K \0 \ 9 + LQ/\d |). From (A4) we 
find 


Xsin(tu 0 t + 4) cos(<o 0 t + <t>)] 



dt a 1 sin(w 0 1 + 4) 


Xcos(tu 0 t + <^)[|sin(fi7 0 ?-(- 0)|] 1 . (A15) 

The integral in the first term of (A15) also occurred inJA6) 
and is zero. The second term may be shown to be zero by 
substituting x = co 0 t + 4 an d splitting the integral, as was 
done in (A7). Thus ^ is a constant for the rigid pendulum 
and may be set to zero by a suitable choice of initial condi¬ 
tions. 

a(t ) obeys Eq. (A3), which gives, with 4 — 0, 

t*2ir/to 0 

a(t) = — K (27r) ~ 1 J dt a 2 col ) sin(<o 0 1 ) | sin 2 (co 0 1 ) 


J ~2.Tr/co 0 * 2 / . \ 

^srnKt) 
o |sin(ft) 0 t)| 

= — 4ATa 2 <i> 0 /(377 r ) — 2 L /(co 0 tt) . 

Equation (A 16) may be written as an integral 

r da (M^ + jk\- i = _r dt , 

Jl = O V (377-) (77“d£> 0 ) ) JO 


(A16) 

(A17) 


which may be solved to give, in conjunction with (A2), an 
expression for 6: 

& — {co 0 )~'{$L /2K ) in cos(w 0 1 ) 

Xtanjtan _1 [a 0 £u 0 (2AT/3Z. ) 1/2 ] 

— 2f7r~ 1 (2iX/3) I/2 ) . (A18) 
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